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Abstract. A summary of the results of a numerical study of the growth of solid particles in the 
vicinity of an azimuthally symmetric density enhancement of a protostellar disk are presented. The 
effects of gas drag and pressure gradients on the rate of growth of dust particles and their settling 
on the midplane of the nebula are also discussed. 



INTRODUCTION 

It has recently been shown that the combined effect of gas drag and pressure gradients 
causes solid objects to rapidly migrate toward the location of maximum pressure in a 
gaseous nebula [1,2]. Such rapid migrations result in the accumulation of objects in the 
vicinity of local pressure enhancement and may increase the rates of their collisions and 
coagulations. In a disk with a mixture of gas and small solid particles, the motion of an 
object is also affected by its interaction with the particulate background. An object may, 
in such an environment, sweep up smaller particles and grow larger during its gas-drag 
induced migration. In this paper I present a summary of a study of the growth-rates of 
dust grains to centimeter- sized objects in a nebula with a local pressure enhancement. 



X : THE MODEL NEBULA 



An isothermal and turbulence-free protostellar disk is considered here with a solar-type 
star at its center. The nebula is assumed to be a mixture of pure molecular hydrogen at 
hydrostatic equilibrium, and small submicron-sized solid particles. The density of the 
gas in this nebula is considered to be 
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In equation (1), M is the mass of the central star, G is the gravitational constant, 
and Kg, T and mjj represent the Boltzmann constant, the gas temperature, and the 
molecular mass of hydrogen, respectively. The quantities po,r m and /3 in equation (1) 
are parameters of the system with constant values. The gas density function as given 
by equation (1) ensures that along the vertical axis, the gravitational attraction of the 
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central star will be balanced by the vertical component of the pressure gradients [2], It 
also has azimuthally symmetric gas-density-enhanced regions on any plane parallel to 
the midplane. On the midplane, the density of the gas maximizes at r = r m . 

The particles of the background material of a nebula are strongly coupled to the gas 
and their motions are only affected by the gas drag. For the submicron-sized particles 
considered here, the rate of the gas drag induced migration is so small that one can 
assume, at any position in the nebula, the gas and particle distribution functions are 
proportional. That is, Pdust(^z) = f Pg{ r ^z), where Pdust^z) is the distribution function 
of the background material. The solid-gas ratio / is taken to be constant and equal to 
0.0034. 



EQUATION OF MOTION 

For an object larger than the submicron particles of the background material, the above- 
mentioned gas-particle coupling is less strong. As a result, the object moves faster 
than the background particles and may collide with them. Such collisions may result 
in adhesion of the background particles to the moving object and increase its mass. 

The rate of change of the momentum of an object due to the sweeping up of the 
background material is proportional to the rate of the collision of the moving object 
with those particles. Assuming the sticking coefficient is equal to unity, one can write, 

dP d£ 

— = n Pd Ur,z)a 2 -(Y-V), (2) 

where a is the radius of the object and d£= \V — U\dt. Quantities V and U in equation (2) 
represent the velocities of the object and the particles of the medium. In writing equation 
(2), it has been assumed that, while sweeping up smaller particles, the density of the 
object remains unchanged and it stays perfectly spherical. It is important to emphasize 
that in this equation, the mass of the object, m, and its radius, a, are functions of time, 

and are related as 

dm 2 d£ 

— = ^p dust (r,z)a-. (3) 



Equation (3) immediately implies 

da _ lf pdast \d£ 
dt ~ 4 V p ) dt" 

where p is the density of the object. 

For small particles such as micron-sized and submicron-sized objects, one can replace 
| V — U| with V re i, the relative velocity of the object with respect to the gas. The radial, 
vertical and tangential components of this velocity are, respectively, given by r, z, and 
r(<p — (Dg), where the motions of gas molecules along the z-axis have been neglected. 
Because the gas is at hydrostatic equilibrium, its angular velocity, 0) g , is slightly different 
from its Keplerian value and is given by 



a _ GM + 1 d£? g (r ,z) 
(r 2 +z 2 ) 3 /2 rp g (r,z) dr 
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In this equation, & g (r,z) = KBTp g (r,z)/m,H is the pressure of the gas. 

An object in the model nebula considered here is subject to the gravitational attraction 
of the central star and the drag force of the gas. The equation of motion of an object in 
this nebula is, therefore, given by 

.. GMm n . , 2 <M 

mR = ~ , ? , ,, ti9 R - ^Pdustl^, Z) a — V re i - F drag , (6) 
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where R(r, z) is the position vector of the particle and 
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represents the drag force of the gas. In equation (7), o is the collisional cross section 
between two hydrogen molecules and A is their mean free path. 



NUMERICAL RESULTS 

The equation of motion of a particle and the growth equation (4) were integrated, 
numerically, for objects with initial radii ranging from 1 to 100 microns. The mass of the 
central star was chosen to be one solar mass, j8 = 1 , r m = 1 AU, and po = 10~ 9 g cm~ 3 . 
The collisional cross section of hydrogen molecules, a, was taken to be 2 x 10~ 15 cm~ 2 , 
and their mean free paths A (cm) = 4 x 10~ 9 /p g (r,z) (gcm~ 3 ). At the beginning of 
integrations, an object was placed at (r,r/10), with a Keplerian circular radial velocity 
and with no motion along the z-axis. Figure 1 shows the growth of two one micron- 
sized objects, one migrating radially inward from (2,0.2) AU, and one migrating radially 
outward from (0.25, 0.025) AU to (1,0) AU, the location of the maximum gas density 
on the midplane. For the comparison, the radial and vertical motions of these objects 
without mass-growth have also been plotted. As shown here, by sweeping up the smaller 
particles of the background, these objects grow to a few centimeters in size and approach 
the midplane in a time much shorter than the time of similar migration without the mass- 
growth. 

Equations (4) and (6) were also integrated for different values of the object's density 
and gas temperature. As expected, the rate of growth of the object increased with an 
increase in the temperature of the gas (Fig. 3). This can be attributed to larger value of 
the pressure gradients at higher temperatures [1]. Increasing the density of the object 
while keeping the temperature of the gas constant had, however, an opposite effect. 
Objects with higher densities tend to grow in size over longer periods of time. The 
inverse proportionality of the rate of the growth of the object to its density can also be 
seen from equation (4). 
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FIGURE 1. Graphs of the radius (left), radial migration (middle), and vertical descent (right) of a 1 
micron object with a density of 2 g/cc. The gas temperature is 300 K. The dashed line indicates migration 
without mass-growth. The middle graphs show an inward radial migration from 2 AU (top), and an 
outward radial migration from 0.25 AU (bottom) to the location of the maximum gas density on the 
midplane, (1,0) AU. The initial value of z was taken to be r/10. 
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FIGURE 2. The growth of a 1 micron object with a density of 2 g/cc for different values of the gas 
temperature is shown on the left. The right graph depicts the growth of similar particle with different 
densities in an isothermal gas with a temperature of 300 K. The particle was initially at (2,0.2) AU. 
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